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^  This  problem  can  be  reduced  to  an  eigenvalue  problem  of  a  fourth  order  ordinary 
differential  equation  in  a  semi-infinite  region.  The  series  that  satisified  the 
boundary  condition  at  infinity  is  asymptotic,  which  cannot  satisfy  the  boundary 
.  condition  at  a  finite  point  and  fails  to  yield  the  eigenvalues.  Theoretica3.1y  this 
difficulty  can  be  overcome,  as  explained  in  this  report,  by  converting  the  divergent 
asymptotic  series  to  a  convergent  factorial  series. 

The  Dartmouth  Time  Sharing  System  (DISS)  was  used  for  the  numerical  computation  of  the 
factorial  series.  It  was  found  that  some  of  the  factorial  series  converge  very 
slowly  and  could  not  yield  any  meaningful  result  even  in  the  utmost  capability  of  the 
DTSS.  It  was  decided  that  the  use  of  the  factorial  series  should  abandoned. 

The  reason  for  the  failure  is  described. 
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1.  Problem  to  be  solved 


We  consider  a  thin  plate  floating  on  -Jpter  stressed  with  uniform 

horizontal  pressure  along  the  periphery  i>f  an  Internal  hole.  We  are 

interested  in  formulating  the  buckling  (trat<iiure  and  the  location  of 

the  maximum  deflection  in  terms  of  the  ':..lus  of  the  hole. 

The  vertical  deflection  of  an  ela;  ■  tc  plate,  which  rests  on  a 

liquid  and  is  subjected  to  a  vertical  v'cd  q  and  the  horizontal 

internal  stress  of  components  N  ,  N  ,  and  N  ,  is  described  by  the 

XX  \ ^  xy 

partial  differential  equation 


DV  w  +  yw  ®  q  +  N 


o  w 

XX  -  2 

^x 


2Jl  “f-  +  N  ^ 
xy  5xdy  yy  ^^2 


(1) 


where  D  is  the  flexural  rigidity,  y  is  the  specific  weight  of  the 
liquid  (ref.  1). 

In  our  problem,  the  deformation  is  cylindrically  symmetric  around 
the  axis  at  the  center  of  the  v.ircular  hole.  Let  r  be  the  radial 


^dr^ 


and  hoop  horizontal  stresses  in  the  plate.  In  the  general  polar 
coordinates 


hole. 

Then  (1)  becomes 

^  +  N  i  H 

[“rr  ^^2  ^  "'ee  r  dr/. 

(2) 

length) 

,  and  N  and  N  ^  are 
’  rr  ee 

the  radial, 

J5: ,  -i-  *wV{«V^-!SJS?*t5WWi?5r:^!r?! 


.2  ,2  ^2 

«  — JT  +  2N  :rr —  +  N  — :r 
XX  xy  f»xdy  yy  ^^2 


rr  gy2  re  hr  \t  hfj  00  V^2  ^^2  r  hr/. 
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V 

r  hp' 


Mansfield  (re£..l)  assumes  that  the  horizontal  stress  N  ,  N  , 

XX  xy 

Nyy  is  in  equilibrium  by  itself  and  shows  that  they  are  derived  from 
a  bih&rmonic  function  cp  by 


» 

“  8y^ 


S  =H 

y'’  Sx^ 


M  =  _  r^iiiiJM.-. 

xy  9x3 y 


In  the  general  polar  coordinates  they  are: 


M  =  i  +  i_  Lift 

rr  r  9r  2  .  2 

r  9p 


N  =  .  ,i  ^ 


re 


9r  r  90.' 


In  our  problem  co  is  a  function  of  r  only  and  must  tend  to  zero 
when  r  becomes  infinite.  Then  we  must  have 


N  =  -  “T  . 
rr  2 

r 


! 


2 


and 


N  „  =  0, 

rp 


P9  ^2  ’ 

where  A  mua*:  be  positive  because  is  pressure.  We  express  them  as 


ay 


N  = 
rr 


and 


N. 


!Lk 

2 


O') 


where  a  is  a  positive  constant  to  be  determined. 
Introudce  the  non-dimensional  length  x, 

^‘viL 


(4) 


Then  (2)  becomes 


4 

,3 

,2 

d  w  ,  2 

d  w 

1-a 

d  w  j  1-a 

— r  +  ~ 
.  A  X 

,  3  ' 

2 

TT +  T" 

dx 

dx 

X 

dx  X 

-  +  w  =  0  (5) 

t 

The  boundary  conditions  we  presently  choose  for  this  problem  are: 
w  =  0 


7~  =  0 

dx 


rat  X  =  X 


o 


w  =  0 


\ 


at  X  w  00 


dw 

dx 


=  0 


where  is  the  value  of  x  at  the  periphery  of  the  hole. 


f 


2.  Asymptotic  Power  Series  Solution 


Wassow  (ref.  2)  describes  a  matrix  theory  of  ordinary  differential 
equations.  According  to  his  theory  the  differential  equation  (5) 
has  an  asymptotic  power  series  solution  of  the  form 


w 


(6) 


where  p^,  X  and  ^  are  constants.  We  will  evaluate  the  constants, 
Wassow’ s  theory,  which  are  concerned  with  proving  the  existence  of 
the  asjTnptotic  power  series  solution,  is  not  convenient  to  this 
purpose.  The  constants  will  be  determined  here  by  d.'rectly 
substituting  series  (6)  into  equation  (5).  Wassow' s  theory,  however, 
is  needed  for  explaining  the  existence  of  the  factorial  series 
solution,  and  will  be  outlined  in  the  next  section. 

The  result  of  the  substitution  of  (6)  into  (5)  can  be  rearranged 


as: 


CO 


n^O 


U,-n 


=  0 


From  A 

o 


0  v/e  find 

4 

X  +1  =  0 


From  Aj^  =  0  we  find 


U  =  - 


For  simplicity  we  put 


1. 


(7) 


4 


Constants  and  P2  are  found,  from  A2  =  0  and  =  0,  as 


p^X  =  -  ^  (a  + 


,2  1  ,  L  ^  9, 

?2^  =  32  (a  ••■  2^  ■*■  2^ » 


respectively.  For  p^  (n  >  3)  we  can  find  from  A^  =  0  (n  >  4)  the 
recurrence  formul: , 


+  (-2(l-a)(n-l)  +  (2n-l)(n-l)(2n-3))  p^_2X  "  + 

+  -^  (2n-l)(2n-5)  (4a  +  (2n-l)  (2n-5)) 


Equation  (10)  enables  us  to  determine  p^X  («  >3)  as  a  polynomial 
of  constant  a  using  (7),  (8),  and  (9)  as  the  initial  conditions.  To 
find  the  polynomial  let 


n  r 

Pn'  ■  “ 


From  the  conditions  (7),  (8),  and  (9)  we  find 


®o  "  ■  8  ’  ®1  “  '  4 


o2__9_  2^5_  2^1. 

Pq  128  ’  ^1  32  »  ^2  32 


1 


Substituting  (11)  into  (10)  we  find  the  recurrence  formulas, 


,  .n  ^n-1 
“4nB  =  6  , 
^n  n-1 


/  *  n  n*  1  c  /  \  A  1  *  o  /  1  \  ^ 

'  "•.,-2  8n-l  ^-2 


.n 


■'"«;-2  -  s'j  ^  «<■»  C\  *  2(.-i>  e“:^  +  g(n)  + 


+  h(n)  (l:l 

-4n0"  =  e“lj  +  f(n)  S"'^  +  2(n-l)  +  g(n)  + 

+  h(n)  0””^  +  h^(n)  0^  ^  (for  n-3  ^  r  ^  1) 

-4n8^  =  f(n)  0^  +  g(n)  +  h  (n) 

where 

f(n)  =  6(n-i)  -  1 

g(n)  =  (n-1)  (4(n-l)^  -3) 
h(n)  =  (n-j)  (n-|) 

For  n  *  3  the  fourth  formula  is  not  needed.  For  n  >  3  the  fourth  formula 

must  be  used  with  r  «  n-3,  n-2  ...,  1. 

4 

Of  the  four  roots  of  X  =  "1»  we  choose 
X  -  (1  ±  i)m 

to  satisfy  the  boundary  condition  of  w  at  x  =  00. 
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The  solution  converging  at  x  ■  e»  is  thiis: 


w  =  e*pf±  i(^  +  f-))  £  s”  a". 


where  C  is  a  complex  constant.  The  real  value  solution  will  be  expressed 


w  =  Aw^  +  Bw2 


with  real  constants  A  and  B,  where 


-x/ir2  r  -n-1/2  /X  .  rfT  \n„n  r 


w 

-x/\f2  r  -n-l/2  /x  ■  nn  \  ^.n 

„2  =  e  jZ.Br 


Or  changing  the  order  of  summations. 


00  GO 


r=o  n-"r 


-K/yfl  ^  r  V  n  -n-1/2  /x  m  \ 

”2  "  Z.  2.  r  j 

r=o  n=r 

The  last  form  of  w^^  and  W2  arc  needed  in  the  next  section  for 
evaluating  constant  a. 
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According  to  Wassow  (ref.  2,  Theor.  12.3)  the  series 

i  +  r)  +  ?))  i  "r 

n=o  r=o 

is  asymptotic.  Still  for  this  case  (15)  is  valid,  because  the  order  of 
summations  in  an  asymptotic  power  series  can  be  exchanged  (ref.  2,  Theor. 
9.5). 

The  asymptotic  power  series  (15),  which  diverges  when  n  is  increased 
indefinitely  for  a  fixed  x,  is  not  convenient  for  determining  constant 
a  and  the  ratio  A:B  from  the  conditions  for  at  x  «  x^.  The  computation 
of  the  asymptotic  power  series  can  be  avoided,  as  explained  below, 
because  the  asymptotic  pw-er  series  can  be  transformed  to  a  convergent 
factorial  series. 

j.  Existence  of  a  Convergent  Factorial  Series  Solution 
A  factoria/.  series  is  a  series  of  the  form 
a 

_ r _ 

x(x+l)  ...  (x+r) 

r=o 

According  to  Wassow  (ref.  2,  Theor.  47.1),  the  differential  equation 
(5)  has  a  convergent  factorial  series  solution.  For  the  proof  he  uses 
a  matrix  expression  of  ordinary  differential  equations.  The  di'. ferential 
equation  (5)  will  be  given  a  matrix  expression  in  this  section  in  order 
to  show  that  our  equation  satisfies  the  requirements  in  Wassow' s 
theorem. 


8 


The  last  four  equations  may  be  lumped  'nto  a  matrix  equation. 


The  entries  of  the  matrix  function  A(x)  are  finite  at  x  =  »  and 
analytic  in  the  neighborhood  of  x  «  a-,  satisfying  a  condition  in 
Wassow's  theorem. 

The  eigenvalues  of  the  matrix  A(oo)  are  the  roots  of  the 
determinant  equation, 


i.e. 

u 

\  +1=0. 

The  eigenvalues  are  therefore  distinct.  Thus  our  equation  satisfies 
all  the  requirements  of  Wassow's  theorem. 

The  resemblance  of  the  matrix  equation  (16)  to  an  equation  of 
single  unknown  seems  to  suggest  a  similar  treatment.  Write 


r: 

t, 


t 

P 


where 


A(x)  =  Aq  +  Aj^x  ^  +  A^x’^, 


^0 


0  0  0  0 

0  0  10 

0  0  0  0 

-10  0  0 

N 


^  = 


"2  = 


0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  -2 

0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  1-a  0 


0  0  0  0 

0  0  0  0 

0  0  0  0 

0  -d-a)  0  0 


The  form  of  A(x)  suggests  to  integrate  (16)  to: 


-1  1  -2 

Y  =  exp  [A^x  +  Aj^  log  x  -  A2X  -  A^x  ] 


(18) 


The  matrix  theory  defines  the  matrix  power  of  e  as 


r=o 


However  the  relation 

A+B  A  B 
e  =  e  e 

is  not  true  unless 

AB  =  BA 

where  A  and  B  are  matrices,  Stcause  of  this  restriction,  (18)  cannot  be 
reduced  to  a  meaningful  form.  The  e  /oidance  of  the  use  ot  (18)  seems  to 
me  the  key  to  the  theory  described  by  Wassow  for  proving  the  existence  of 
the  solution  in  the  form  of  (6). 
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4.  Various  Expressions  of  x  ^ 

The  convergent  factorial  series  solution  of  the  differential 
equation  (5)  can  be  found  by  substituting  the  factorial  series 
expressions  ofx^(p=l,  2,  3,  ...)  into  the  asymptotic  power  series 
solution  (12).  The  existence  theorem  (ref.  2,  Theor.  47.1)  guarantees 
the  validity  of  the  procedure. 

In  this  section,  following  the  traditional  treatment,  various  forms 
of  factorial  series  expressing  x  ^  (p  =  1,  2,  ...)  will  be  derived. 

In  the  numerical  computation  all  forms  of  the  factorial  series  developed 
here  were  tested  but  none  could  yield  satisfactory  results  in  all  the 
computations. 

The  proof  starts  from  recalling  the  definition  of  the  gamma 
function, 


e“^  d?, 


defined  for  p  >  0.  Changing  the  variable  §  to  t  by  ^  =  xt,  we  find 


-P  1  r”  -xt  ^p-1 

'  =  Fw  J  "  ' 


Changing  the  variable  t  to  z  by  z  =  e  ,  we  have 

iiiii::!  r 

r(p)  J, 


,-p  =  ■  r  jjc-i  ^^p-i 


(19) 


Let 


cp(z) 


("1/ _ 

r(p) 


(log  z) 


p-1 


(20) 
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If  we  substitute  the  Taylor  series  of  cp(z)  developed  around  z  =  1, 

to 

CP(Z)  ^  (£>^'’•'(1)  (z-1)”, 

n=o 

into  (19)  we  find 

00  ^ 

n=o 


where  B(x,  n+1)  is  the  beta  function 

.1 

B(x,  n+1)  =  I  z^“  (I-z)”  dz  (22) 

*^0 

Various  factorial  series  expressions  of  x  ^  stem  from  (21)  as  shown 
below. 

A  sophisticated  device  is  necessary  to  derive  the  successive 
derivatives  (p^”^(l):  Let  G(log  z)  be  an  analytic  function  of  log  z. 

One  sees  easily  that 

n-1 

^  G(log  z)  =  z‘"  Y  (-1)''  (log  z)  (23) 

dz  ^ 

v=o 

for  n  >  1  with  certain  constants  s"  which  are  the  same  for  all  functions 

=  V 

G(5) ,  where 


G^*"^  (log  z)  =  G(r)l 

••dl  '€  =  log  z 

••  cF  tl 

Let  G(f)  =  e  'to  determine  S  ,  where  c  is  a  Cynstant.  Then  G(log  z) 

=  Z  ^  and  d'^/dz'^  G(log  z)  =  (-1)"  c(c+l) . . .  (c+n-l)z  ^  G^^^  (log  z)  = 


13 


Hence  (23)  reduces  to 


n-1 


c(c+l) . . . (c+n-1) 


-L 


v=o 


For  n  =  1  and  2,  (24)  easily  yields 

s\  =  1,  s]  =  0 


2  2  2 
4  =  ^1  =  ° 


For  a  large  n,  recurrence  formulas  for  computing  S  will  be  found 
Increasing  n  in  (24)  by  1  yields 


n 

c(c+I) . . . (c+n)  =  ^ 

v=o 


gn+1  ^n+l-v 
V 


I 


Comparing  the  two  formulas  we  find  the  relation, 

n  n-1 

„n+l  n+l-v  /  ,  \  ^  oil  V 

S  c  =  (c+n)  >  Sc 

V  4*  V 

v=o  v=o 

Because  c  is  arbitrary,  the  following  must  hold; 
=  1 


,-.n+l  „n  .  .n  _  ,  , 

S  =  S  +  n  S  , ,  for  1  <  v  <  n-1 

V  V"!  =  = 


V  V 


and 


=  nS”  ,  =  nJ 

n  n-1 


n+1 


(24) 
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15 


v^^  -*r .  •  "'/T v^-v '■* /yv.^* 


,-P 


=  Y  ^  S”  B(x,  n+l), 

^  n;  n+l-p  ’ 

n=p“l 


(26) 


Substituting 


B(x,  n+l)  = 


rCx),  nl 
r(x+n+l) 


(26)  yields 


,-P 


I 


,n 


n+l-p 


x(x+l) . . .  (x+n) 


(27) 


n=p-l 


For  p  =  1  the  right-hand  side  reduces  to  1/x  because  =  1  a’  =  0 
for  n  >  1.  Therefore  (27)  is  effective  for  p  >  2. 

Let  0  be  a  positive  number  and  transform  (22)  as  follows: 


.1 


B(x,  n+l)  =  [*  ^  z  (1-z)” 


dz 


=  f  (l-(l-z))"®  (1-z)”  dz 

^  f  ^  (-1)®  f’'^)  (1-z)®’^"  dz 


s=o 


=I  rnf 

imm  n 


x+p“l  /I  vS+n  , 
z  (1-z)  dz 


s=o 


L/d+s-1  (n+s) ! _ 

'  s  /  (x+D)...(x+p+r+s) 


s=o 


Substituting  this  into  (26),  we  find 
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'  V, 


n=p-l 


(x-fp ) . . .  (x+D+n) 


where 


r=p-l 

For  p  =  1,  (28)  simplifies  considerably.  In  the  following  we 
derive  the  same  result  by  a  simpler  method.  We  now  start  from  the 
formula 


1  r  x-1  . 

7=1  2  dz 


which  transforms 


r  x+p-l  -p  , 

=  z  z  ^  dz 


J  (l-d-z))""  dz 


r  x+p-i  ,,  .s  , 
z  (l-z)  dz 


=  ^  B(x-tp,  s+1) 


/i:>+s-l'\  _ 

'  s  /  (x+p).. 


(x4p+s) 


Thus  we  find 


i  ^  , _ E _ .  _ eifiii} - + 

X  x+p  (x+p )  (x4p+l)  (x'hj)  (x+o+1)  (x+p+2) 


*/ 
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By  a  similar  method  a  multiplier  on  x  can  be  introduced  into  the 


denorainctors  of  any  of  the  formulas  (27),  (28),  and  (29).  This  is  the 


form  used  by  Wassow  (ref.  2,  Theor.  47.1).  We  however  do  not  need  such 


a  generalization,  and  will  not  derive  such  a  formula.  In  the  following 


we  rather  derive  formulas  in  which  x  is  raised  into  a  numerator. 


Transform  (22)  as  follows: 


B(x,  n+1)  =  f  (l"(l-2))*“^  (l“z)”  dz 

4 


=  f  (-1)''  f  (1-r.)”'*'’'  dz 


r=o 


,  JL..  .  r 

n-t-1  ^  rJ  (n+r+1) 

r=l 


Substituting  this  into  (26),  we  find 


k'P  -  bP  -I-  7  bP  ,  (30) 

o  t  v: 

r=l 


where 

w  „Vi 

5?  =  r 

o  L  (n+DJ 

n^p"! 


Vy?  =  V  — 

r  n! (n+r+1) 

n-p"l 
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For  p  =  1,  (30)  simplifies  considerably.  In  the  following  we 
derive  the  same  result  by  a  simpler  method.  We  start  from  the  integral 


Which  transforms 


i=r  (l-Cl-z))"""^  dz 


=  y  (-1)’'  J  (1-z)’'  dz 


•I 


Thus  we  find 


X  “  ^  ^  2!  31  4!  •“ 


or,  writing  similary  to  (30), 


X  o  ^  r  r! 


where 


,  1  1 
^r  r+1  ’ 


r  >  0 


The  right-hand  sides  of  (30)  and  (31)  are  a  form  called  by  Norlund 
(ref.  5)  as  a  Newtonian  interpolation  series 


:-x)  U-x) . . .  (r-x) 
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The  range  of  convergence  of  a  Newtonian  interpolation  series  is 
quite  different  from  a  power  series;  if  the  former  converges  at  x^,  it 
converges,  similarly  to  a  factorial  series,  for  x  >  x^.  The  strange 
impression  one  may  have  on  reading  this  statement  will  reduce  when  one 
observes  that  (32)  ends  with  a  finite  number  of  terms  for  any  positive 
integer  x. 

The  convergent  coordinate  X  is  the  least  x^  .•:uch  that  for  x  <  x^  - 
X  the  series  does  not  converge.  It  is  given  by 

n-1 

( 

log 


_  lim  sup 


n-> 

when  X  >  0  and  by 


_  lim  sup 
n  oo 


s=o 


(33) 


log  n 


Vlog  n 


(34) 


s=n 


when  X  <  0.  The  notation  lim  sup  c  means  that  if  c  is  less 

k  — ..A  «  ’^■’■1 

h-r  oO 

than  the  preceding  numbers  c^,  c^,  ,  the  number  should  be 

discarded  in  taking  the  limit. 

The  range  of  convergence  of  a  factorial  series 
00 


a  s] 


x(x+l) . . . (x+s) 


(35) 


s=o 


is  similar  to  (32).  Its  convergent  coordinate  X  is  given  by 


l™  Clog  I  r  a  )/U.g  : 

n-^co  ®  ' 


(36) 


s=o 
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when  y  a  diverges  and  by 


s=c 


)/log 


s=n+l 


when 


^  a^  converges.  In  (36)  X  ^  0,and  in  (37)  X  <  0.  The 


s=n+l 


formulas  of  convergent  coordinates  in  the  above  are  listed  from 

n 

Norland  (ref.  4  and  5).  Borel  (ref.  3)  proves  that  the  Newtonian 
interpolation  foriiula  converges  or  diverges  when  the  correspond inpr 
factorial  series  converges  or  diverges.  (Note;  For  clarification  of 
the  last  statement  I  must  re-oheck  Borel  (ref.  3),  which  was  borrowed 
from  M.I.T.).  This  theorem  is  important  for  substituting  a  Newtonian 
interpolation  series  for  a  factorial  series. 

Using  (25)  we  can  show  that 

r 

(l-x)...(r-x) 

v=0 

Use  of  this  formula  reduces  (30)  to 


09 


n“0 

where  we  can  find  that 
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-p+1 


cP=(.i)^r  -1 .  X  (r\ 

n  ^  ^  L  (^+1)1  L  r-n  7-r-p+l 

^n+p-1  r=n 

For  p  =  1,  the  coeficients  reduce  to 

c;  =  (.u»£ 


(e+i) :  7-n 

-n+p- i 

Eq  (39)  is  the  formula  to  be  used  to  fird  a  power  series  expansion  at 

X  =  0  as  the  analytical  continuation  of  the  asymptotic  power  series 

solution.  Coefficients  c^, 

0 


“  I 


■P+1 


=P-1 


(t-rU  7-r-p+l 


r»o 


diverge,  although  the  divergency  is  not  obvious  except  c^, 


€=0 


If  the  analytical  continuszion  to  x  =  0  of  an  asymptotic  power  series  is 
convergent,  we  must  be  able  to  rearrange  the  order  of  summation  to  find 
a  finite  value  at  x  =  0.  (i"u:  Viumerical  computation  was  not  successful, 

however,  in  this  oasr  because  of  the  slow  convergence.)  Formulas  (30)  and 
(31)  are  interesting  in  itself  because  it  shows  close  relationship 
between  a  factorial  series  and  a  Newtonian  interpolation  series.  Moreover 
(38)  is  useful  for  finding  the  successive  derivatives  of  (1-x) . . . (r-x) . 
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5.  Determination  of  Constant  a  and  the  Ratio  A;B 

Constant  a  and  the  ratio  A:B  in  the  as^nnptotic  power  series  solution 
are  found  from  the  loundary  conditions  at  -  x^.  The  asymptotic  power 
series  diverges  when  n  ia  increased  indefinitc-ly  for  a  fixed  x.  It 
must  be  converted  to  a  factorial  series  or  a  Newtons  an  Interpolation 
series  whose  convergences  arc  guaranteed.  Such  a  sei Les  can  oe  found 
in  our  problem  by  substituting  a  suitable  expression  of  x  ^  in  the 
previous  section  into  the  asymptotic  power  series  solution  and  the  order 
of  summatiops  is  changed  accordingly,  oubstitution,  however,  mrv  be 
postponed  until  algebraic  manipulations  for  solving  the  simultaneous 
equations  for  constant  a  and  the  ratio  A:B  are  finished  because  the 
algebraic  manipu  .ations  are  admissible  to  an  asymptotic  power  series. 

From  (12)  we  find 


g  -  ,-3  ^  )1^  v;  a'. 


n 


n=o 


r=o 


where 


Y°  =  -1,  =  0  for  r  >  1 

'o  'r  = 


n  „n 


Yn  = 


-e. 


Yr  = 


=  -b"  +  (7  '  n)  b”"^  for  0  r  r  <  n-1 


In  real  form 


,  dw,  dw„ 

dw  ^  .  __1  ..  R  _1 

dx  dx  dx  ’ 


(40) 
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I 


T.th 


dWj 

dx 


-n- 


r  ^  “  2 


n-o 


‘“172  I  \ 

r=o 


4. 

4 


^"'2  ^-x/)f2  'C  ^"-2  ^  .^x  ^  inJJn. 


dx 


=  -e 


n=o 


n 


)t ''" 


r=o 


Or  exchanging  the  order  of  summations, 


(41) 


dw 

dx 


00  00 


1  -x/\r2r‘  rr-  n 

-=-e 

r=o  n=r 


•n-2 


cost  ^  +  iHiiir-  ) 
‘^°®Sr2  4  J 


(42) 


dw 

dx 


00  00 


2  -xm 

=  -e 


Zrr  n  """2  .  /X  .  (n-])n\ 

a  Yr  ^ 


r=o  n=r 


Use  of  (13)  and  (40)  in  the  condition  at  x  =  x^  yields 


B  w, (x  )  V  dx  /  / ^  dx  / 

1  o  X  /  X 

oj  o 


(43) 


Eliminating  the  ratio  A:B,  the  equation  for  determination  of  constant  a 
is  f  ’and 


t=o 


(44) 


where 
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wl  « 


-p  lP 


(45) 


p=t 


and 


r+p-t 


I 


n  .p-n  .  (p“2n+l)n 

V  P  sin 
'r  t-r  4 


(46) 


r=o  n=r 

When  p=t  constants  lP  approach  zero  as  p=t  increases,  and  for  a 
fixed  t  constant  increases  rapidly  in  their  absolute  values  as  p 
increases  indefinitely.  The  wide  range  of  values  covered  by  constants 
lP  was  one  of  the  causes  of  the  computational  difficulty. 


Using  (27),  functions  Fj.(Xq)  expressed  as  convergent  factorial 


series, 


and 


(47) 


for  t  >  2.  The  convergence  of  these  factorial  series  is  guaranteed  because 
w  and  dw/dx  expressed  as  factorial  series  converge. 
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When  approaches  zero  as  the  limit,  we  replace  with 


whose 

values  G^(o)  at  x  =  0  for 

n  >  0  are  finite: 

00 

r+1 

G„(0)  = 

[  ir 

•  I  ^0^ 

r=l 

P=2 

G^(0)  =  lJ  +' 

00 

r+1 

■  1  "-f  ^r-p+l 

(48) 

r=l 

p=2 

00 

r+1 

G,(0, 

1  ' 

[  ^P.p« 

r=t-l 

p=t 

Therefore  constant 

a  tends  to  a  certain  value  when 

X  approaches  zero  as 

0 


the  limit. 


When  x^  increases  indefinitely  we  may  use  the  asymptotic  power  series 
(45)  for  evaluating  constant  a  from  (44).  We  can  show  in  this  case  that 


k  =  a/x  (50) 

o 


approaches  a  certain  constant  as  the  limit  when  x^  increases  indefinitely, 
where  k  is  the  root  of 


^  lJJ  k"  =  0  (51) 

k=o 

It  is  noted  that  the  asymptotic  power  series  (45)  cannot  be  used  at  a 


finite  x^  for  evaluating  constant  a. 

The  series  (48)  converged  rapidly  and  the  computation  of  constant 


a  in  the  neighborhood  o*"  x^  =  0  was  carried  out.  The  result  is  shown  in 
Table  1.  The  series  (51),  however,  did  not  numerically  converge,  because 
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i 
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I 

i 

t 


f 
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the  number  k  was  so  large  that  the  series  could  not  converge  even  in  the 
utmost  capability  cf  the  DTSS. 

The  most  serious  limitation  of  the  DTSS  was  the  narrow  range 

(approximately  from  lE-39  to  lE+38)  of  the  magnitude  of  the  numbers 

00 

that  it  can  handle.  The  graph  ^  ^n  touches  the  x  axis 

k=o 


extremely  gently  and  the  root  k  of  equation  (51)  could  not  be  easily 
discovered. 

In  order  to  determine  A/B  in  (43)  computation  of  the  middle  term. 


//Win 

-W2(x^)/w^(x^).  and  the  right  term,  -(--j  /  . 

X  /  X 


was  carried  out 


Of  o 

through  25  terms.  It  was  found  out  that  the  two  series  did  not  converge 
sufficiently  for  all  values  of  x^  to  test  whether  or  not  they  were  equal 
(Table  2). 

Further  computations  were  not  attempted. 

The  language  used  on  the  DTSS  was  the  BASIC. 


i 
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Table  1.  Values  of  a  as  function  of  x  . 

0 


X 

0 

a 

k  =  a/x 

0 

3.18364 

lE-6 

3.18364 

lE-5 

3.18368 

lE-4 

3.18405 

lE-3 

3.18775 

IE'2 

3.22414 

lE-1 

11.1805 

1 

17.1648 

17.16 

2 

48.9837 

24.49 

3 

68.3702 

22.79 

4 

84.172 

21.04 

5 

99.2279 

19.84 

6 

113.586 

18.93 

7 

127.65 

18.24 

8 

141.59 

17.69 

9 

155.48 

17.28 

10 

fails  to  converge 

